Cosmological phase space analysis of the F(X) — V(<f>) scalar field and bouncing 

solutions 
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We analyze the dynamical system defined by a universe filled with a barotropic fluid plus a scalar 
field with modified kinetic term of the form L = F(X) — V ((/>)■ After a suitable choice of variables 
that allows us to study the phase space of the system we obtain the critical points and their stability. 
We find that they reduce to the ones defined for the canonical case when F(X) = X. We also study 
the field energy conditions to have a nonsingular bounce. 



I. INTRODUCTION 

Scalar fields play an important role in cosmological 
models because, due to their simplicity and adaptabil- 
ity, they can account for different interesting phenomena. 
They are some of the most popular choices for modeling 
cosmological scenarios such as inflation [l[ and dark en- 
ergy [2(, and they also have been studied in the context 
of dark matter models [H, bounce cosmology [3 and dif- 
ferent unification models of those phenomena @ . 

The proposal that the Lagrangian could be a general 
function of the kinetic term X — —g IJ,l 'df 1 4'd v 6j2 and 
the field <fi was introduced in cosmology first in [6J in the 
context of inflation and then used for dark energy models 
in Q . Different particular forms of the Lagrangian C — 
p(X, </>) has been studied for different reasons [2| . 

In this paper we study the general class of models with 
sum-separable Lagrangian C — F(X) — V((f>). Several as- 
pects of this type of scalar fields have been studied in the 
literature. Phenomenology of the inflation models arising 
from them p| M, topological defects [l(| , supersymme- 
try extension [ff|, boson stars [13, unification models 
of dark matter and dark energy [5|, [f3| and unification 
of dark energy, dark matter and inflation [Til . [To] ] . This 
model also offers the possibility of being understood as 
a vacuum energy density V coupled with a barotropic 
fluid |f6[. It reduces to the canonical scalar field when 
F(X) = X. 

In section|TI]we study the system of autonomous differ- 
ential equations related to this class of scalar fields. This 
method i s eq uivalent to the one used for canonical scalar 
fields in [T7| and allows us to identify the general be- 
havior of the cosmological solutions associated with the 
present Lagrangian. This method has been applied to a 
wide range of cosmological models, for example [l8l - l26j ]. 
It can be used to determine the presence and stability of 



solutions of cosmological interest, such as those with de 
Sitter phases or with scaling behaviors. 

One possible application for this type of Lagrangian 
is the generation in the early Universe of a nonsingular 
bounce, in which the state of the Universe goes from 
collapsing to expanding for a^O. The bouncing models 
have been proposed as alternatives to inflation [53, [H| 
and as a way to evade a singular big-bang, as in the pre - 
big-bang scenario [1^] , in the ekpyrotic universe (30l - [33j j , 
or in other multifield models |34| . 

In section ULTl we study how, when the density of this 
field is allowed to be negative, it can drive a bounce. As 
the variables defined for previous the dynamical system 
analysis are not suitable to study this phenomenon, we 
redefine the system as in Ref. [34| in order to study this 
case and obtain the conditions to accomplish a bouncing 
behavior. To obtain the bounce, the scalar field has to 
violate the Null Energy Condition (NEC) possibly giving 
rise to instabilities. Some works have been made trying 
to erase these instabilities with ghost condensate scalar 
fields [HI, HI HI (however see Ref. H3). Here we will 
only consider the dynamics of homogeneous cosmologies. 



II. AUTONOMOUS SYSTEM FOR 

C = F(X) - V(</>) 

For a spatially flat Friedmann-Lemaitre-Robertson- 
Walkcr (FLRW) cosmology filled with a scalar field with 
a Lagrangian of the form £ = F(X) — V(4>), where 



x = --d^d^, 



(i) 



and a matter component with density p rn and equation 
of state p m = u m p m , the equations of motion are 



1 



[2XF X -F + V + Pm ] 



(2) 



* josue@ciencias.unam.mx 
t jorge.cervantes@inin.gob.mx 
t david.wands@port.ac.uk 



H 



[2XF x + {\ + w m ) Pm ], 



(3) 



2 



where H is the Hubble factor. These equations can be 
combined to imply the conservation of the total energy 
momentum tensor. We can suppose additionally the con- 
servation of the scalar field and barotropic fluid energy 
momentum tensors separately, which is the case when 
there is no interchange of energy between the two com- 
ponents. In that case the barotropic component satisfies 



the equation p m oc a 



-3(l+i* 



for a constant u> m , where 



a is the scale factor, and the scalar field satisfies 



dN 



(2XF X -F + V) + 6XF X = 0, 



(4) 



where the subindex X means differentiation with respect 
to that variable. The time differentiation here has been 
changed to dN — d log a, a variable that for an expanding 
FLRW model can be used as the independent variable 
instead of the cosmological time, with the relation dN = 
Hdt. 

In order to obtain the autonomous system we define 
the variables 



x 



where x 2 is proportional to the kinetic part of the energy 
density 



p k = 2XF X - F . 



(6) 



and y 2 to the potential part of the energy density py — 
V. They are equivalent to the ones used in the analysis 
for canonical scalar fields pjj ■ We will also need to define 
the auxiliary variables 



F 

2XF X - F '■ 



2X 



3\2XF X - F 



■sign(>), 



(7) 
(8) 



where the former corresponds to the change in time of 
the potential, as can be seen if we write it as 



M P1 dlogV 



V 3 \pk\ 



dt 



(9) 



and the latter corresponds to the equation of state for 
the kinetic part of the Lagrangian, as the kinetic part of 
the pressure is P k — F. In the case of a canonical scalar 
field F(X) = X and the auxiliary variables turn out to 
be uj k = I and a = ^/2/3\ for V oc e - A< ^ Mpl as defined 
inRef. [13 • 

The equation of state of the scalar field can be obtained 
in terms of the new variables as 



W0 



y 



x z + y z 



(10) 



The evolution equations for the first two variables (x 
and y) can be written as 



dx 
dN 

dy 
dN 



2XF XX + F 



A" 



dX 



H 



2^/3(2XF x - F)M Pl H dN H 2 



H 



2VWMpiH 2 ' H 2 



(11) 



(12) 



The common H / H 2 factor can be obtained from equation 
([3]) dividing by H 2 and replacing the original for the new 
variables 

H 3 

where we made use of the equation ^ in the new vari- 
ables 



y 



1 ■ 



(14) 



Now in order to calculate the first term in the evolution 
equation (fl"2"|) we only have to substitute the values of the 
new variables 



2VSVM Pl H 2 



-^oxy. 



(15) 



For the first term in equation (|lll) we use the continuity 
equation (j4]) that can be written as 



dX 
dN 



3F 



{2XFxx + Fx)u k 



ay_ 

x 



(16) 



so that the evolution equations ([TTl[T2l become, in terms 
of the new variables, 



= — \oy 2 — xiuju + 1)1 + 

dN 2 1 y 



dy 
dN 



-x [(1 +u) m ){\ - y 2 ) + x 2 (uj k - uj m )] 
3 

-~ayx + 
3 

-y [(1 + u> m )(l - y 2 ) + x 2 (uj k - oj m )] 



(17) 



(18) 



The evolution equations for the variables uj k and a can 
be obtained using the equation (I16[) and the definition of 
X. But we have to define new auxiliary variables that 
depend on the second-order derivatives of the Lagrangian 
potentials. The evolution equations are 



~dN 



1 2Euj k +uj k -l 
2E + 1 



u k 



X 



da 
dN 



-3a 2 x(T - 1) + 

3a(2E(u k + 1) + oj k - 1) i a ,, 



(19) 



(20) 



2(2S+l)(w fc + l) 



LO k + 1 

x . 



3 



where the auxiliary variables are defined as 
„ XFxx 



r 



Fx 

vv^ 
vl 



(21) 
(22) 



The new second-order derivative variables T, S will 
have evolution equations in terms of the dynamical vari- 
ables and new third-order derivative variables, and so on. 
In order to truncate this succession of equations we can 
consider fixing the functions F(X) and V (</>). 

The first of these assumptions is to choose the potential 
related variable T as a constant. For it to happen we need 



for T 7^ 1 or 



V^)=V ^-4>o) ini - T) 



(23) 



(24) 



for r = 1 . The second assumption is to consider the case 
in which 



F{X) = AX 71 , 



(25) 



where A and ij are constants, in this case u>k = 1 / (2?y — 1) 
and equation (fP9)) is trivially satisfied. In the following 
we will use these assumptions. 

The dynamical system will be reduced to Eq. (fl"T| for 
the evolution of x, Eq. (fT8|) for the evolution of y, and 
an equation for the evolution of a that due to the choice 
of F as a power law becomes 



da , , 



3cr(l - Ul k ) 

2(l+w fe ) 



x 



(26) 



A. Critical points 



The autonomous system of equations written above 
can be analyzed if we consider its critical points, in which 
the equations ([XT]) and (fI5]l are equal to zero, correspond- 
ing to x and y constant. In the first instance we will not 
consider the evolution equation (|26|) . but if the critical 
points (xo,yo) depend on a they will not be truly fixed 
unless we ensure that a is constant. 

The variables x 2 and y 2 correspond to the fraction of 
the energy density contained in kinetic and potential en- 
ergy of the scalar field, as can be seen from (TT4"1) . The 
condition of constancy for the critical points implies that 
these variables have a constant contribution to the total 
energy density, which can happen in three scenarios: (i) 
if x 2 + y 2 is equal to one, meaning that all the energy 
density comes from the scalar field, (ii) if they are zero, 
meaning no contribution, or (iii) if they are between zero 
and one, corresponding to what is also known as a scal- 
ing solution, meaning that the energy density of the field 



scales at the same rate as that of matter. The three be- 
haviors are of cosmological interest and are present for 
the critical points of generally defined parameters x and 
y if they satisfy equation (fl"4"|) . 

To present the critical points, we have labeled with 
Latin letters those that reduce in the canonical case to 
the ones studied in [38(, and with Greek letters to the 
ones with no correspondence. 

For x a = and y a = this corresponds to the scalar 
field not contributing to the energy density of the uni- 
verse. 

If x = and y ^ the equation of state becomes — 
— 1 that is an interesting case from the cosmological point 
of view due to the possibility to describe dark energy or 
inflation phenomena. In this case the evolution equation 
for x reduces to 



dx 
dN 



2 ay 



(27) 



which requires a = 0, a potential that doesn't change 
with time. On the other hand the evolution equation for 
y reduces to 



(28) 



that can be zero for cj m = — 1 or y = 1. Both cases cor- 
respond to a FLRW model filled with fluid with equation 
of state —1. 

• In the first = and y a = 1 , the Friedmann 
equation in the new variables ([T4"l) implies that the 
matter field has zero energy density and the only 
component of the model is the <fi field. 

• In the second case xp = and yp is arbitrary, then 
there are contributions from the barotropic fluid as 
well as from the scalar field. 

If y = and x ^ the potential energy is zero and the 
equation of state reduces to = Wk- The energy density 
of the field is stored in the kinetic part. The evolution 
equation for y vanishes and the one for x reduces to 



dx 3 . 

= — Xiuju 

dN 2 V k 



■u) m ){x - 1), 
that can become zero in two cases: 



(29) 



• For Xb = 1, yb = this corresponds to the density 
of the model coming entirely from the kinetic part 
of the field 4>. 

• For uj m — uj^, with y 7 = and x y arbitrary, the 
equation of state of the field is the same as the 
equation of state of the matter. It corresponds to 
a kinetically driven scaling solution. This type of 
solutions are important in cosmology because in the 
case of dark energy they have been proposed to 
alleviate the coincidence problem This case, 
however is not completely what in the literature is 



4 



called a scaling solution in the sense that it can 
only reproduce a constant equation of state of the 
matter when the Lagrangian of the field satisfies 



F(X) = AX {1+ ^' 2 ^ 



(30) 



For example if the energy density of the matter 
satisfies a relativistic equation of state, we need 
F(X)= AX 2 such that u k = 1/3. The latter hap- 
pens in the unified dark matter models based in 
Scherrer's Lagrangian F(X) = F Q + F m (X - X ) 2 . 
It is known [39( that for high energies in which 
X ^> Xq the model can have a radiation-like be- 
havior and this is because F and X can be disre- 
garded, approximating to (pO)) . 



The last case is when both x and y are different from 
zero. From (Ti~7l [18]) we can see that the critical points 
satisfy 

X= h( U3k + 1± + l) 2 - W) ■ (31) 

In this case there are two different critical points, the first 
one has the form 



u k + 1 ' 



Vc 



y/(cj k + l) 2 - a 2 
Wfc + 1 



(32) 



It corresponds to a cosmology filled with the scalar field, 
as can be seen from (|14[) which in this case corresponds 



to x 2 + y 2 = 1 with zero matter density. The equation of 



state of the system will be 



W0 



1 +LO k 



- 1 



(33) 



The second nonzero critical point corresponds to 



Xd 



lid 



1 



\J (uj m + l)(wfc - bJ m ) 



(34) 



where the equation of state in this case is uj^ — u) m , in 
other words corresponding to a scaling solution. In the 
canonical case with exponential potential we will recover 
the scaling solution of Ref. [ljj ■ The fraction of the total 
energy density stored in the scalar field will be 



2 _ (1 + 0Jk)(l + Urn) 



(35) 



The points defined in (|32|) and (134)) depend explicitly 
on <7, that in general is an evolving quantity. It means 
that, unless the variable a is also fixed, those points won't 
be critical points of the system. Setting the evolution 
equation for a equal to zero gives us the condition T = To 
with 



2(1 + w fc ) 



(36) 



This relation between the derivatives of potential and 
kinetic terms in the Lagrangian has to be accomplished 
in order to have the critical points (13"2l IMf . For the 
canonical case, as u>k = 1 then the critical points are 
fixed only for r = 1, that from (j2~4"]l corresponds to the 
exponential potential, as expected. For the noncanonical 
case as Wfc ^ 1 then there will be a relation between the 
exponent in the kinetic and the one in the potential term 
of the form C = AX r > - B{<j> - <fo) n with 



>1 



(37) 



if this relation is not satisfied, the critical points won't 
be truly fixed. In the appendix we show that when the 
system satisfies this relation, it is invariant under a set of 
symmetry transformations, which turn allows to reduce 
the number of degrees of freedom. The same symme- 
try invariance happens for the canonical scalar field with 
exponential potential as proved in Ref. [4lj . 

The stability of the critical points can be analyzed by 
the matrix of the derivatives of the right hand side of 
equations (TTTl IT5|) . Analyzing the eigenvalues of the ma- 
trix we obtain the results of table |U 

The critical point (a) presents a behavior of unstable 
node for u)k < u) m which can drive the scalar field density 
towards bigger values even if it starts with small density. 
The saddle point behavior that was already obtained in 
the canonical case is recovered here when Lu m < u>k ■ Point 
(a) corresponds to slow roll behavior as a = and the 
potential dominates, and it can be a saddle point or a 
stable node depending on the equation of state of the ki- 
netic part. Point (b) can be stable, unstable or a saddle 
point. In the canonical case the stable behavior is not 
obtained. Points (c) and (d) have the same stability be- 
havior as in the canonical case except that the conditions 
get modified by u>k as stated in the table. The lines (/3) 
and (7) are obtained when the equation of state of matter 
is the same as that of the kinetic part or the potential 
part of the Lagrangian and can be stable or unstable. 
The cosmological relevance of these solutions is further 
discussed in the conclusions. 



It is interesting to point that, except for the canonical 
case, the Lagrangians studied here with uik and T con- 
stants can't be reduced to the case C — Xg(Xe x ^), that 
is considered in Ref. as the general form for a scalar 
field with scaling solutions. The difference from the case 
studied there is that we are not considering a coupling 
between the field and the barotropic fluid as in their case. 



B. Critical points at infinity 

The former critical points a, b, c, a, /3 and 7 corre- 
spond to the situation in which the dynamical variables 
are finite. This is always the case for x and y as Eq. (TT4"|) 
requires both variables to be smaller or equal than one. 
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TABLE I. Stability and existence of the critical points assuming —1 < u> m < 1- The points labeled with Latin letters reduce 
in the canonical case to the ones already studied in the literature [38l |. the points with Greek letters are new. 





X 


y 


Existence 


Stability 






(a) 








Always 


Unstable node for cj^ < lu 7JI 
Saddle point for u m < uj^ 







(a) 





1 


a = 


Saddle point for oJk < —1 
Stable node for > — 1 


i 


-1 


(b) 


1 





Always 


Unstable node for > {^ m , cr — 1} 
Stable node for uJk < {^mi cr — 1} 
Otherwise saddle point 


i 




(c) 


IT 


j{l + u k ) 2 -* 2 


r - r 


Saddle point for a 2 > (1 + Wfe)(l + Wm) 
Otherwise stable node 


i 


a 2 

1 




1 + Wfe 


V (l+"fc) 2 


<r(l + w fc ) > 
fr 2 > (l+o; fc ) 2 




(d) 


1 + U) m 


S/(UJ k — UJ m ){l + U)m) 


r = r 


cr — ci;^ + yixJ m j 
Stable node for < 1 


(1 + UJ k )(l + U>,„ ) 




(7 


(7 


U>rn < 

a > v/(l + Wfc)(l+Wm) 


8(l + w fc )(l + w m ) 2 
Stable spiral otherwise 


<7 2 


(P) 





Arbitrary 


cr — and — —1 


Stable line for UJk > — 1 
Unstable otherwise 


2 

V 


-1 


(7) 


Arbitrary 





w fc — w m 


Stable line for xa > 1 + UJk 
Unstable otherwise 


x 2 





However for a we can see from the definition (JT)) that it 
can become infinity as pk the kinetic energy density or V 
the potential tend to zero. To study this case we make 
a change of the variable a to E = 1/a and study the 
possibility of it becoming zero. 

Considering uik and T constants, the evolution equa- 
tions CEZ1 [H [26]) become 



and the factors 



dx 3 ,y 2 , 



-x [(1 + w m )(l - y 2 ) + x 2 (uj k - w m )] 



dy 
dN 



dE 
dN 



3xy 
' 2E 



(38) 



(39) 



-y [(1 + ui m ){l - y 2 ) + x 2 (w fe - w m )] , 

Mr- D- ^1(^ + 1)=-^ | (40) 



In order to have a critical point at E = it's required 
that the terms y 2 /E and xy/E each vanish. These factors 
can be computed considering that as cjk and T are con- 
stants the kinetic term is a power- law F = AX n , and the 
potential term is either a power-law V — B<j) n or an ex- 
ponential V — Ce~ A */ Mpi . For the power-law potential 
the variable E has the expression 



nMpi 



^"^ sign^A^'- 1 )/ 2 , (41) 



ir_ 

s 

E 



n5 



MIpiH 2 




3M Pi H 2 V 3(2r?- 1)A 



sign 



-sign 



n-2)/2j^l/2 



n-l x (l- n )/^ 2 ^ 



(43) 



In order to have these two terms equal to zero at the 
same time as E — ¥ 0, we requiere = and n > 2. With 
these conditions the variable y becomes zero too and the 
evolution equations for x and y get reduced to 



= — x(x 2 — 1)(ujl 

dN 2 V A 

dy 



dN 



= 0, 



(44) 
(45) 



which implies that the critical points occur when x — 
0, x — 1 or UJk = w m . These three cases correspond 
to the already studied critical points (a), (b) and (7), 
respectively. 

The second case happens when the potential is an 
expone ntial V — Be~ Ac ^/ Mpl . In this case however 
V = \AVPc cannot become zero, except for the trivial 
case in which B = 0. This means that j/ 2 /E diverges as 
E — » 0, which implies that (f39| cannot be zero and there 
is no critical point at infinity. 



C. Analysis of the phase space 

In this subsection we plot the phase space defined by 
the equations (fTT]) . (fT8|). (|26|) for several potentials. An 
important case happens when the potential satisfies equa- 
tion (|36[) . For example this occurs for a Lagrangian of 
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y 




0.0 ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 X 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 1. Projection of the phase space along the (x, y) plane 
for a Lagrangian of the form £ = AX 2 + B/4> 4 , with initial 
condition a = 1.5. The solutions tend to the critical points (c) 
and (d) studied in table [I] As a changes, the critical points lie 
on the light green line for the (c) and the red segment of circle 
for (d). This behaviour can be better seen in the figures [5] and 
[3] corresponding to different projections of the same system. 
The critical point curves are plotted with the same colors. 

the form C — AX 2 + B / (j) 4 . In this case the extra param- 
eters have values of uik — 1/3 and r = 5/4. In Figs. [TJ [2j 
and [3] we plotted the two dimensional projections of this 
system for us m — for an initial condition of a — 1.5. We 
can see how the solutions approach to the critical points 
(c) and (d) from table Q] depending on the values of a. 
The critical points lie on a curve in the 3-dimensional 
phase space and due to the evolution of a the solutions 
tend to different points in those curves. 

As we have stated, for the canonical scalar field the 
condition (|36[) implies an exponential potential. In that 
case a is a constant determined by the exponent in the 
potential, as V cx e -V3^/(V2M Pl ) < Tlie pnase spaC e in 
that case is effectively two dimensional. In Fig. |4]we can 
see the behaviour for a = 1.5 with oj m = 0. In this case 
the solutions approach only to one point, as there is not 
evolution in a. 

When the condition (pIB")) is not satisfied, we don't have 
the critical points (c) and (d), but we can still plot the 
system. For example, for the Lagrangian C = AX 2 +Bcf> 2 
in which case the values of the auxiliary parameters are 
uj k = 1/3 and V = 1/2. In Fig. [5] we plotted the three 
dimensional phase space system for solutions that start 
with a — 1.5. We can see that the system evolves towards 
big values of cr, this happens because a oc l/(f> in this case, 
and the system goes towards small values of 4>- I n fact, 
it crosses <f> = in a finite time, in which the variable a 
is not useful to describe the system. We can see that the 



(J 




X 



0.0 02 0.4 0.6 0.8 1.0 



FIG. 2. Projection of the phase space along the (x,a) plane 
for a Lagrangian of the form £ = AX 2 + B/<f> 4 . The system is 
the same as in figures [1] and [3] The solutions where chosen to 
start in a = 1.5 and for a constant x they evolve in different 
directions due to the different values in y. See the explanation 
in Fig. [1] 

solutions don't tend to any critical point. 



III. BOUNCE COSMOLOGY 

In this section we consider a nonsingular bounce (a ^ 
0) in a FLRW cosmology filled with the scalar field 
C = F(X) — V(<fr) and a barotropic fluid with constant 
equation of state u) m . For a bounce to happen we need 
the evolution of the scale factor to go from decreasing to 
increasing as a function of time. In terms of the deriva- 
tive of the scale factor this implies that at the bounce it 
has to satisfy a(tf,) = and d(tb) > 0. The first condi- 
tion can be translated in terms of the Hubble parameter 
as Hb = 0, but this means that the dynamical variables 
defined in equations ([5]) in the last section will diverge. 
Besides that, the independent variable N that we have 
used to parametrize the evolution of the system is no 
longer well defined at the bounce, as d/dN — H~ 1 d/dt. 
Those complications arise from the fact that our choice 
of dynamical variables was adjusted to study a cosmol- 
ogy with increasing a. Accordingly, in order to study a 
bouncing FLRW metric we need to define new variables 
adapted to the current problem. 

Also we have to note that the total energy density 
of the model is zero at the bounce, which can be seen 
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log(x) 



FIG. 3. Projection of the phase space along the (x, a) plane 
for a Lagrangian of the form C = AX 2 + B/cj> 4 . The system is 
the same as in figures [T] and [2] The solutions where chosen to 
start in a = 1.5 and they evolve towards different directions 
depending on the values of x and y. See the explanation in 

Fig. in 



y, 




FIG. 4. Phase space for the canonical scalar field with expo- 
nential potential and a = 1.5. The phase space for this system 
is two dimensional, unlike in the case of nonexponential po- 
tentials or noncanonical kinetic terms with three dimensional 
phase spaces. The solutions tend to the critical point (d), 
scaling solution. 




FIG. 5. Phase space for the Lagrangian £ = AX 2 + B<p 2 . The 
solutions don't tend to any critical point. They evolve towards 
high values of a because this parameter is proportional to <j)~ 
for this Lagrangian. 



from the Friedmann equation ([2]) . If we suppose that the 
energy density of the barotropic component is positive 
p oc a~ 3 ( 1+u,m > then the energy density of the field has 
to be negative, something that will be considered in the 
definition of the dynamical variables below. 
Now let us define the new set of variables 



V 



sign(V) . 



(46) 



where p k is given in Eq. ([6]) and the absolute values come 
from the fact that we are interested in the behavior of 
both, positive and negative energy densities. The new 
independent variable defined in analogy to N is 



dN 




3Mp\ 



(47) 



The evolution equations for the above variables can be 
then written as 

dx 3 

-7T T = -oK^fc - W m )sign(p fc ) + (1 +UJ m )(x 2 - y\y\)] 
dN £ 
3 

+-x [(uj k + 1)5 - cry\y\sign(p k )] , 
dy 3 

-TT T = T.V [~<7 + (wjfc + l)Z - ay\y\sign{p k )] . (48) 

In general we also need to evolve a, its evolution equation 
can be obtained from equation (|20[) , transforming to the 
new variables as 



da 
dN 



-3cr 2 (r- 1) + 

3a(2E(uj k + 1) +uj k - 1) 
2(2H + l)(w fc + l) 



(49) 



((ojk + i)^ - vy 2 ) , 
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FIG. 6. Schematic projection of the phase space for a non- 
canonical nonphantom system with pk > and uik = — 5 and 
a ~ a/2/3. The bounce occurs when the solutions cross the 
vertical (thick, red) line. 



FIG. 7. Schematic projection of the phase space for a field 
with pk < 0, LJk = 1/6, a ~ —a/2/3 and w m = 1/3. The 
bounce occurs when the solutions cross the vertical (thick, 
red) line. There is no purely kinetic (pi/ = 0) bounce. 



The above variables are well behaved only for pk / 0, so 
neither of the possible cases of purely potential bounce, 
nor a change in sign for pk after the bounce will be stud- 
ied here. 



A. Conditions for a bounce 

In the following we work with the phase space defined 
by the set of equations l|48p . Due to the relation between 
the dynamical variables (x, y) and (x, y) from the previ- 
ous section, the critical points of both systems coincide. 
It's easy to check that the points of table Hare also crit- 
ical points of the new system with the transformation 
x = 1/x and y — y/x, except for those with x — in 
which the new variables diverge. 

In this subsection we will use equations (|38"]) to study 
the evolution of the systems near the bounce. In gen- 
eral we have to consider also equation (|5P|) to make a 
representation of the three dimensional phase space as 
in the Fig. [5] of the previous section, however we will 
not consider this equation because we are interested in 
the behaviour only close to the bounce and the variable 
a won't evolve much during this short time. For this 
reason the plots of the phase spaces of Fig. |H] and [7] 
that will be studied with more detail in this section, cor- 
respond only to schematic representations of the phase 
space near the bounce. For figures [8] and [9] the represen- 
tation corresponds to the actual phase space because for 
those Lagrangians a is a constant. 

Besides equations the system has to satisfy the 



Friedmann constraint (|2"]h which translates into 

* 2 - V\v\ - O m = 1 X sign(p fe ), (50) 

where VL m = p m /\pk\ corresponds to a dimensionless den- 
sity parameter for the barotropic fluid component. As 
Q m is assumed to be nonnegative, we obtain the expres- 
sion 

x 2 - y\y\ > 1 x sign(p fc ), (51) 

which defines the allowed regions of the phase space. For 
the nonphantom case pk > 0, the Friedmann constraint 
becomes 

x 2 - sign(I/)y 2 > 1, (52) 

which for y positive corresponds to the region inside the 
branches of the hyperbola x 2 — y 2 = 1, and for y negative 
to the region outside the circle defined by x 2 + y 2 = 1 , as 
can be seen in figure |U In the pk < case the condition 
(jST) translates into 

x 2 -sign(V r )y 2 > -1, (53) 

which for y > corresponds to the region below the 
hyperbola y 2 — x 2 — 1. For y < this condition is 
satisfied for all the values, as we can see in the figure [8] 
From the definitions (|46| we can see that 

• x > corresponds to the regime of an expanding 
cosmology, 

• x < corresponds to a contracting cosmology, 
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X 



FIG. 8. Phase space (x,y) for the special case of a phantom 
system F(X) = -X and V oc e' X(p/Mpl (such that cr is con- 
stant) plus a barotropic radiation component u) m = 1/3. The 
bounce occurs when the solutions cross the vertical (thick, 
red) line. The spiral in the graph corresponds to the critical 
point (d) studied in the previous section. 



• x = corresponds either to a bounce, a recollapse 
or a static cosmology. 

For the case of x = we can use the information con- 
tained in the derivative to study whether we are dealing 
with a bounce or a recollapse: 



dx 
dN 



> corresponds to a bounce, 



• < corresponds to a recollapse, 

• = gives not enough information and one has 
dN . ° . 

to consider higher derivatives or analyze the neigh- 
boring phase space. 

To see which of the above cases occurs in the phase 
space of our system we use x — in the evolution equa- 
tions (J3HJ) - I n particular for the evolution of x we obtain 



dx 
dN 



[(w fc - u m ) sign(p fe ) - (1 + Lu m )y\y\] . (54) 



As we stated above this expression has to be positive for 
a bounce, which implies a condition in the parameter y 



V > 



OJk — u ri 



1 + UJ„ 



sign(p fe (wA; - uj m )) . (55) 



In addition we also have the condition ([Si]) for the case 
x = 



y 




2 10 1 X 



FIG. 9. Phase space for the special case of a canonical scalar 
field with potential V oc e~ x ^^ Mpl (such that a is constant) 
plus a barotropic radiation component ui m = 1/3. All the 
solutions that cross x — have negative dx/dN which corre- 
sponds to recollapse. The bounce is not possible. 



y < -1 X sign(p fe ). (56) 

To analyze the above conditions we hrst suppose > 
0. In this case the inequality (|55| transforms to y < — 1 
and (|55|) to y 2 < (co m — Wfe)/(1 + w TO ). For the two 
conditions to be satisfied in an interval of y is necessary 
to have Uk < —1- For example, in the case of a canonical 
scalar field, one has F(X) = X and consequently pk > 0, 
but as uik = 1 the system of a barotropic component 
and a canonical scalar field cannot give rise to a bounce, 
as was already shown in (42J. This can be seen in the 
phase space of figure [9] in which all the solutions that 
cross the y axis move from positive to negative values of 
x, corresponding to recollapse. 

For pk < the conditions for the bounce become 



\ — — sign(w m - uj k ) < y < 1 , (57) 

y i + w m 

which can be satisfied for an interval of y as long as 
ujk > —1. The original phantom field with F(X) = —X 
satisfies OJk = 1 and, as shown in figure [51 can have a 
bounce behavior. 

In order to have a purely kinetic bounce, in other words 
one with y = the conditions above state that the den- 
sity of the scalar field pk has to be negative and ujk > uj m . 
Figure [7] shows a case in which the later is not accom- 
plished and then there is no purely kinetic bounce. 
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Bounce 



FIG. 10. The densities of the barotropic fluid (blue, dash- 
doted line), the scalar field (red, doted line), and the total 
density of the universe (green, continuous line) respectively 
as a function of the scale factor. The total energy density 
tends to zero at the bounce and for smaller values of a is 
negative, which is forbidden. 



The two conditions in the previous paragraph can be 
generalized. First, to obtain a bounce one needs the to- 
tal energy density of the field to be negative in order to 
compensate for the positive barotropic energy density in 
the Friedmann equation 



H 2 = 



1 



3-Mf?! 



(P4>- 



= 0, 



(58) 



where p^ = py + pk is the total energy density in the 
field. Moreover, the total equation of state of the field 
W0 has to be bigger than that of the barotropic fluid in 
order to have a positive energy density for a > atouncei 
as can be seen in figure [TO] Otherwise we will be dealing 
with a system that exhibits positive energy density only 
for a < cn, ounce corresponding to a recollapse. 

The above two conditions are in fact the same as 
those in the expressions ([S3)) and (151)1) in terms of the 
dynamical variables. For the first one, the negativ- 
ity of the energy density pk + py can be translated as 
1 x sign(pfc) + p v /\pk\ < or from the definitions of the 
variables ([irJ|) as 



y < -1 x sign(p fc ), 



(59) 



which corresponds to expression (|56[) . The condition on 
the total equation of state of the field, in terms of the 
dynamical variables can be written as 



Uk - y|y|sign(p fc ) 
1 + y|y|sign(p fe ) 



(60) 



which can be transformed into ()55[) after some algebra 
and using the expression (|5T))) . 

The conditions on the field to have a negative energy 
density and an equation of state greater than that of 
matter implies a violation of the Null Energy Condition 




FIG. 11. The bottom left region corresponds in the p — p 
plane to the part which can drive a bounce, with p < and 
o)0 > oj m with uj m = 1/3. The upper right region is the 
one that satisfies the Null Energy Condition. The dash-doted 
lines cannot be crossed by k-essence Lagrangians like the ones 
consedered here 14311. 



(NEC) that states that p0 + be positive, as can be 
seen in Fig. 111! In the last years an extensive literature 
has been produced studying fields that violate the NEC. 
The main reason for that interest is because the cur- 
rent measurements of the dark energy equation of state 
slightly favor models with oj^e < —1 0' However fields 
violating the NEC might have several types of instabili- 
ties, for example imaginary sound speed which results in 
an increase of inhomogeneities in small periods of time 
[43|, 53, or decay of the vacuum into negative energy par- 
ticles of the field plus positive energy particles |45T - [47I |. 
The inclusion of higher order terms in the Lagrangian 
has been proposed as a method to obtain particles with 
positive energy in the so called Ghost condensate models 
[3TI l35l [H, |4§] , however usually these extra terms add 
new stability problems to the models, and is not clear if 
there is a well behaved high energy theory to account for 
them . Due to those problems, a recent series of works 
has been published studying fields in which the introduc- 
tion of certain symmetries can ensure the stability of the 
model in spite of breaking the NEC (50l - [5^ |. However 
those models, so called Galileons, have dynamics which 
was not studied in this paper. 



IV. CONCLUSIONS 

As we have seen, the system of equations ^ [3]) for the 
Lagrangian C = F(X) — V((f>) can be rewritten in terms 
of the dynamical variables as (|171 IT5|) . This system 
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allows us to understand the dynamical behavior of the 
universe under different initial conditions. The critical 
points and their stability are summarized in table [H This 
system is naturally adapted to study Lagrangians with 
kinetic terms of the type F(X) oc X v and potentials 
V{<fi) oc i^ 1 ^ 1 "" 1 ) or V{4>) oc e _A< ^ such as those studied 
for k-inflation in The canonical case and its critical 
points are recovered for F(X) — X. 

In general the critical points (a), (/3), (7), (c) and (d) 
are present only for particular choices of the Lagrangian 
C = F(X) — V ((/)), as happens for the canonical scalar 
field in which the points (c) and (d) are only present for 
exponential potentials. The conditions for their existence 
are summarized in table HI 

The point (a) corresponds to the slow roll scenario in 
which the potential dominates (y = 1) and its derivative 
is zero (a = 0). The case with ujk < — 1 is interesting for 
inflationary models as it corresponds to a saddle point, 
offering an explanation of how the universe could enter 
in the slow roll regime and exit eventually. For that case 
it is also necessary to study the dynamical behavior of 
a to understand the conditions for it to evolve towards 
zero, something that was not analyzed in this paper. 

The potential dominated line (/?) has the cosmologi- 
cally interesting behavior of an equation of state of —1, 
however it requires that the barotropic fluid has the same 
behavior, something that is very restrictive. 

The kinetic dominated line (7) corresponds to critical 
points of the system only when ujk = w m , for example if 
F(X) oc X 2 when the barotropic fluid is radiation. It can 
happen for example in the purely kinetic unified model 
studied in Ref . [391 ] , in which the proposed Lagrangian be- 
haves as a radiation fluid for high energies. An interesting 
extension to this purely kinetic model is the addition of 
a potential term to the Lagrangian which could leave the 
kinetic dominated line stable at early times, setting the 
initial conditions necessary for a later evolution as dark 
matter plus dark energy if the potential becomes flat at 
late times, see also Ref. 151]. 

The scalar field dominated solution (c) and the scal- 
ing solution (d) are not in general critical points of the 
system except for the case in which the potentials in the 
Lagrangian satisfy the particular relation (|36|) . This re- 
lation for the canonical scalar field means that the po- 
tential has to be exponential, and for the noncanonical 
field means that the Lagrangian has to satisfy ([3"T]) . As 
in the canonical case, the scaling solution corresponds to 
a stable node or a stable spiral, and the scalar field dom- 
inated solution behaves as a stable node or saddle point. 
If condition f|3~B]) is not satisfied, even if dx/dN = and 
dy/dN — for a particular time, the variables will evolve 
because the time dependence on a will drive x' 7^ and 
y' 7^ as time passes. However in the cases in which 
is satisfied we can obtain scaling solutions despite 
the fact that the Lagrangian cannot be reduced to the 
form C = Xg{Xe x ^), which is studied in Ref. gjj as 
the general form of scalar fields with scaling solutions; 
but we are not considering here an interaction with the 



matter component as in that case. 

In order to study a bouncing cosmology we had to 
redefine the dynamical variables to some more suited to 
the problem as (|46[) . We obtained the conditions (|55l [56]) 
necessary for a bounce. In the phase space it is seen as 
the possibility to have a crossing of the y axis from the 
negative to the positive x region. 

The dynamical variables (x, y) and (x, y) are related by 
the transformation x = 1/x and y = y/x which means 
that the critical points of both dynamical systems coin- 
cide when both are valid. This happens when pk and pv 
are positive, otherwise the variables x, y are not defined, 
and when i^O. It can be seen that the points of table U 
are also critical points of the new system except for those 
with x = 0. 

We split the analysis of the bouncing system in two 
cases, pk negative (phantom scalar field) and pk positive, 
and obtained that in order to have a bounce we need 
ujk > — 1 for the first case and Uk < —1 for the second 
one. 

For a canonical scalar field, we know that a negative 
potential can lead to a crossing of the x axis (H = 0) 
only for recollapse, and not for a bounce. Here we showed 
that for certain values of uik a bounce is possible even for 
Pk positive, giving the possibility of a potentially driven 
bounce. We also showed that the conditions ([551 [55)1 
obtained in terms of the dynamical variables can be ul- 
timately understood as p^, < and uj^ > tj m , better 
seen from figure [TO] as the conditions to have zero en- 
ergy density at the bounce and positive energy density 
immediately after and immediately before it. 

We showed that the field has to violate the Null En- 
ergy Condition in order to account for the bounce, as can 
be seen in figure [IT] This is a well known result that can 
have implications concerning the stability of the field. It 
this paper we didn't deal with the inhomogeneous per- 
turbations, however it has been argued that this type of 
Lagrangians have both classical and quantum stability 
problems when they violate the NEC [13, IH- All the 
former arguments make us conclude that possibly fields 
as simple as F-V are not good candidates to violate NEC 
and therefore to produce a bounce. The study of other 
types of fields might be in order but it escapes the pur- 
pose of the present paper where only the homogeneous 
dynamics of the fields was considered. 
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Appendix: Symmetry for particular Lagrangians 

The critical points (b) and (d) from table U exist only 
for canonical scalar fields with exponential potential or 
for scalar fields whose Lagrangians are of the form 



with 



C = AX r > - B( 



V 



(A.l) 



(A.2) 



In these cases the system presents a symmetry that al- 
lows the number of degrees of freedom to be reduced to 
two, and the dynamical system to be described only by 
x and y. For the canonical scalar field with exponential 
potential this symmetry was described in 41]. 

The equations of motion plus the continuity equa- 
tion for the barotropic component can be written for a 
Lagrangian of the form (|A.1[) as 



H 2 = 



1 



n d JL = 

dN 

dprn 



[(27? - 1)A¥" + Bcj> n + p m 
[2 V AX^ + {l+u m ) Pn 



2M 2 Y 



dN 
d 



-3(l + u m )p v 



dN 



((2?7 - 1)AX^ + Bc/) n ) = -&r]AX^ , 



(A.3) 

(A.4) 

(A.5) 
(A.6) 



where for simplicity we considered tf>a = 0. Here <f>, X, 
and p m are the independent variables and the transfor- 
mation 



X -> i 2n X . 

Prn * s P? 



(A.7) 



will leave invariant the equations of motion as long as the 
Hubble parameter also transforms as H — > £ nTI H, but its 
transformation is already determined by the relation 



which implies that H transforms as (; n ~ 2ri H . In order to 
have the correct transformation relation for the Hubble 
parameter then it is needed that nr\ — n — 1r\ which is 
equivalent to the relation (|A.2I) . only in that case the 
transformation (|A.7|) will represent a symmetry of the 
system leaving invariant the equations of motion. 



The presence of the symmetry transformation (|A.7j) 
when (|A.2j) holds means that the number of degrees of 
freedom in the equations of motion can be reduced by 
one. For this a set of variables invariant under the trans- 
formation needs to be defined, in this case x and y are 
already invariant. Any dynamical variable can be writ- 
ten in terms of those two variables, for example a satisfies 
the relation 



x 

s I - 



2/r, 



(A.9) 



where s is a constant defined by the parameters in the 
Lagrangian as 



, = - xi -M Pl nB 1 ^ n (A(2r 1 - l))" 1 / 2 " . (A.10) 



From this relation, the dynamical system can be rewrit- 
ten as 



dx 3 r /y\2/(«*+i) 

dN = 2 [sXyi 
3 



x(tu k + 1)] + (A.ll) 



[(1 + uj m )(l - y 2 ) + x 2 (u> k - LJ m )] , 



dy 
dN 



(!) 



- — SX 

2 \x 

-y [(1 + w m )(l - y 2 ) + x 2 (u> k - oj m )] 



(A.12) 



X 



2 ( H dN 



(A., 



corresponding to only two equations for two variables. 
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